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Abstract 


This  report  summarizes  a  preliminary  result  contained  in  an 
upcoming  paper^  documenting  continuing  research  effort  in  bubbly 
ocean  surface  scattering.  As  such,  details  will  be  referred  to 
that  paper.  This  report  is  written  for  a  timely  comparison  with 
different  sets  of  experimental  data  and  is  hoped  to  shed  some 
light  on  a  possible  reason  why  there  were  apparent  differences 
among  experimental  data  sets. 

Since  the  ONR  Special  Research  Program  (SRP)  has  publicized 

the  importance  of  bubbly  ocean  surface  effect  on  acoustic 

scattering,  many  papers  dealing  with  the  subject  have  appeared. 

The  model  described  in  this  paper  is  based  in  part  on  a  NATO 
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presentation  delivered  in  May  1992. 
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1.  Introduction  to  the  Model 


Consider  an  incident  plane  wave  velocity  potential  4^  - 
exp[  i  {  *x-(»jt )  ]  entering  the  wind  generated  bubbly  surface  layer 

from  below  (see  Fig.  1),  where  -  k(a,(3,-y),  k  -  incident 
acoustic  wave  number,  (a,(3,-y)  -  incident  wave  direction  cosines, 
w  -  circular  acoustic  frequency,  t  -  time,  x  «  (x2D,z)  -  (x,y,z) 

■  3D  position  vector,  x2D  ■  (x, y)  -  horizontal  position  vector, 
and  z  »  mean  depth.  The  objective  is  to  estimate  the  total  field 
of  reflected  waves  coming  down  from  the  bubbly  layer.  Due  to  a 
rough  ocean  boundary  and  fluctuating  sound  velocity  of  the  layer 
medium,  there  will  be  both  surface  and  volumetric  scattering. 
Within  the  single  scattering  theory,  there  are  two  chains  of 
alternating  volume  and  surface  scattering  events. 

The  first  chain  starts  from  the  volumetrical  scattering  of 

the  incident  wave.  Then  the  up-going  volumetrically  scattered 

waves  will  be  surface  scattered,  and  volumetrically  scattered 

again.  As  described  in  reference  3,  single  surface  scattering 

can  be  symbolically  characterized  by  a  multiplication  factor  of 

(1+dA+dAdA*)  where  dA  is  a  Fourier  Stieljes  component  of  the 

surface  roughness  h/  i.e.,  h(x2D^  “  I  exP  t-ik2D*x2D]dA(k2d) . 

The  single  volumetric  scattering  can  be  shown  to  be  symbolically 

characterized  by  a  multiplication  factor  dr,  which  is  a  Fourier 

Stieljes  component  of  (Sc/Cq),  the  fractional  acoustic  velocity 

2 

fluctuation  i.e.  2k  (Sc/Cq)  -  Jexp( -ik *x )dT( k )  - 

Jexp( -ik • x ) T( k )dk .  Because  both  dA  and  dr  are  considered  small 

perturbation  parameters,  the  scattered  field  generated  by  the 


alternating  repeated  scattering  events  will  eventually  become 
negligible  -  beyond  second  order  in  dA  and  dr. 

The  second  chain  starts  from  a  rough  surface  scattered  wave 
♦  .  According  to  the  single  volumetric  scattering  described  by 

O 

4 

Batchelor  ,  the  incident  wave  will  penetrate  the  bubbly  layer 
approximately  unchanged  despite  the  ongoing  volumetric  scattering 
de '•cribed  above.  The  scattered  waves  ♦  will  be  volumetr ically 
scattered  and  followed  by  the  similar  alternating  scattering 
processes  as  in  the  first  chain  until  the  scattered  waves  become 
negligible.  The  total  reflected  and  scattered  field  is  the 
combination  of  all  the  down-going  waves  generated  by  every 
surface  and  volumetric  scattering  events  or  processes  described 
above.  The  detail  of  obtaining  the  total  reflected  and  scattered 
field  can  be  found  in  reference  1.  In  this  report,  major  steps 
taken  in  reaching  a  solution  will  be  summarized. 

For  each  volumetric  scattering,  the  first  Born's 

4 

approximation  as  presented  by  Batchelor  is  assumed  appropriate  - 
the  single  volumetric  scattering  approximation.  Previous 
investigators  McDonald^  and  Henyey®  utilized  the  same 

4 

approximation.  In  addition,  they  also  utilized,  as  Batchelor 
did,  the  far-field  approximation  for  an  individual  plume. 

Approach  here,  however,  do  not  identify  individual  plumes.  A 
continuous  layer  (without  an  abrupt  mean  impedance  discontinuity) 
of  an  infinite  horizontal  extent  is  assumed  to  exhibit  random 
sound  velocity  fluctuations  induced  by  those  plumes  which  appear 
randomly  in  space  and  time.  Volumetrically  scattered  waves  are 
divided  into  two  groups.  The  first  group  propagates  away  from 


the  sea  surface  and  contributes  directly  to  the  reflected  field. 
The  second  group  propagates  toward  rough  sea  surface  and  surface 
scattered.  Since  volumetric  scatterers  are  near  ocean  surface, 
these  up-going  volumet r ically  scattered  waves  should  be  in 
near-field  forms. 

For  each  sea  surface  scattering,  the  reflected/scattered 
field  is  obtained  by  applying  a  perturbation  method^  to  a 
pressure  release  boundary  condition.  Most  of  the  reflected 
energy  contributes  directly  to  the  reflected  field  below  the 
bubbly  sea  surface  layer.  A  small  part  will  be  volumetr ically 
scattered  as  described  above. 

The  total  reflected  field  4>r  ( x2D,  z ,  t )  consists  of  all  above 
down-going  (away  from  ocean  surface)  waves.  To  the  second  order 
in  dr  and  dA,  there  are  a  total  of  nine  component  waves.  These 
component  waves  are  re-grouped  according  to  the  perturbation 
order.  It  can  be  shown  that  four  of  these  component  waves 
without  the  surface  roughness  combine  into  the  basic  results 
before  far-field  approximation  utilized  by  McDonald5  and  Henyey5 . 
At  this  stage,  the  total  reflected  field  is  expressed  in  terms  of 
random  variables  dr  and  dA  as  a  result  of  turbulent  wind  action. 
Therefore,  the  total  reflected  field  takes  different  values  from 
one  realization  to  another  and  only  statistical  averages  are 
relevant.  An  appropriate  statistical  analysis  is  described  in 
the  following  section. 


/ 


2 .  Statistical  Analysis  and  the  Backscattering  Strength 

Predictive  Model 

Before  applying  an  appropriate  statistical  analysis, 
relationships  between  (Sc/Cq)  and  other  measurable  physical 
parameters  should  be  investigated.  This  approach  will  make  it 
possible  not  only  to  evaluate  consistencies  of  different  physical 
parameter  measurements,  but  also  to  enable  a  use  of  particular 
physical  parameters  that  are  easier  to  measure  for  scattering 
strength  predictions.  As  shown  in  reference  1,  these 
relationships  can  be  summarized  into  one  general  form: 

( Sc/c0 )  — Fx  exp(-pz)X(x2D)  (2-1) 

...  in  which  exponential  depth  decay  with  decay  constant  3  is 
assumed.  X  represents  one  of  random  surface  characteristics  in  a 
horizontal  space  such  as  fractional  acoustic  velocity 

fluctuation  (5c/c0)s,  number  of  bubbles  per  unit  volume  nq, 
roughness  element  height  fraction  of  sea  surface  covered  by 

white  caps  W,  and  acoustic  scattering  cross-section  per  unit 
volume  Mvq.  is  the  corresponding  coefficient  for  each  x 

selected.  Expressions  of  these  coefficients  are  given  in  the 
following.  The  minus  sign  in  equation  (2-1)  came  from  an 
expression  relating  (Sc/Cq)  to  N  at  a  depth,  see  Clay  and 
Medwin^.  Though  more  complicated  expressions  can  be  similarly 
derived  for  a  given  bubble  size  spectrum,  following  expressions 
assume  an  existence  of  a  dominant  bubble  size  of  radius  a. 


wR  -  resonant  circular  acoustic  frequency 

g  -  gravity  acceleration 

aCH  "  -  Charnock  constant 

T  -  surface  tension 

/j  m  dynamic  viscosity  of  water 

C  -  minimum  surface  wave  velocity 

m  ■* 


This  means  that  the  volumetric  scattering  characteristics 
dT( k )«r( k )dk  can  be  estimated  not  only  from  the  direct 
measurement  of  (8c/Cq)s  but  also  from  the  indirect  measurements 
of  Nq,  W,  Mvq-  The  proportional  relationship  between  (Sc/cQ) 

and  requires  the  assumption  of  Charnock' s  dimensional  power 

O 

law  to  be  3  which  is  higher  than  his  original  value  of  2  and 
also  higher  than  the  maximum  empirical  value  of  2.5  suggested  by 

O 

Wu  .  Although  there  is  no  direct  evidence  that  the  assumption  is 
justified  at  this  time,  it  is  a  plausible  way  by  which  the 


roughness  element  roughness  and  surface  roughness  b  can  be 
simply  correlated.  Depending  on  the  assumption  of  controlling 
physical  parameters  in  Charnock's  relationship,  F^k  takes 
different  expression.  F  ^  results  if  gravity  only  is  important 
(Charnock's  original  assumption).  f^2  results  if  gravity, 
surface  tension,  and  viscosity  are  all  important  (suggested  by 

O 

Wu  ) .  results  if  gravity  and  surface  tension  are  important 

as  suggested  by  reference  1.  In  the  following  preliminary 
estimates,  F  ^  is  utilized.  The  reason  for  using  this  expression 
as  the  first  trial  in  checking  experimental  data  is  that  gravity 
wave  spectrum  is  the  only  one  documented  at  this  time. 

The  simple  depth-wise  exponential  decay  in  equation  (2-1)  was 

confirmed  by  measurements  of  different  physical  quantities, 
g 

Thorpe  had  shown  also  that  the  functional  form  is  a  solution  to 
a  vertical  diffusion  equation.  The  spatial  coordinates 
dependence  of  (Sc/Cq)  in  a  product  form  of  equation  (2-1)  implies 
also  the  following  product  form; 

r(k)  «  r(k2D,kz)  «  rz(k2)  rx(k2D)  (2-2) 


where 

2nrz(kz)  -  -2k2  Jexp(ikz  zQ-fizQ ) Fx( zQ )dzQ 


(2-3) 


(2n)2rx(k2D) 


Jexp(  iji2D— 2D }  2D }  d— 2D 


(2-4) 


The  extent  of  F^  dependence  on  depth  depends  on  a  selected 
physical  variable  for  X-  One  interesting  observation  made  by 
Wu10  was  the  fact  that  the  bubble  size  spectra  of  two  experiments 
were  found  to  be  approximately  invariant  with  either  the  depth  or 
the  wind  velocity.  Therefore  depth  dependence  of  a  is  probably 
negligible.  There  is  a  small  depth  dependence  for  in  the 

following  estimates,  however,  the  value  of  F  will  be  assumed  a 

A 

depth  averaged  constant.  Thus,  integration  in  equation  (2-3)  can 
be  easily  performed.  r  and  Fw  are  considered  determinestic 

z  x 

while  rx  and  x  are  considered  random. 

Now,  the  reflected  field  velocity  potential  $r ( x2D, z , t )  can 
be  re-written  in  terms  of  random  variables  dr  and  dA.  An 

A. 

appropriate  statistical  average  is  the  cross-correlation 

★ 

<+r(x2D'°'t>  ♦r<x2D*i.°.‘>> 

where  <  >  represents  an  ensemble  average  while  £  represents  a 

horizontal  separation.  The  reflected  field  spectrum  n(k_r  2d)  is 
then  given  by  the  Fourier  integral 

n(kr,2D)  "  (V2"*2  JexP(  ikj. ,  2d*-^}  <^*r(— 2D'0)^r’lr(— 2D  +  L'0)>dL- 


The  result  is  given  below. 


n(kr,2D}  "  t  Specular  Terms]  S(k.f2D  -  kr,2D> 


+  (i/4M2*rJ(I(ks-krB)/k„»2«r8l(k^-k„)/k„^xx(ki  2D-krf2D) 


+  4  kiz*AA  (-i,2D~-r,2D) 


+  ki2[2nrzI<kZ“krz)-/krz+2nrzl(kz-krz)/krz^xA^i,2D-kr,2D^ 


(2-5) 


where  kr  -  ( iir  f  2D'  krz ) '  -i  “  (ki,2D,kiz>'  -  *  (k2D'  kz 


). 


r«I<kz> 


-r  <-k.  -k  )  +  r  (-k.  +k  )  +  r  (k.  -k  )  -  r  <k.  +k>, 


z'  " iz  z’  ~z'  "iz  "z'  'z'  "iz  z'  z'  "iz  ”z 

and  specular  terms  -  1-  energy  loss  due  to  scattering  in  all 

directions.  Spectra  and  *XA 

correlation,  X  correlation,  and  x-h  cross-correlation, 

respectively  are  defined  below. 


#  for  roughness  (h) 
AX 


<h(x2D)n*(X2D+l2D>> 

<X(x2D)X*(X2D+^2D)> 

20^  2D+— 2D^ 

Jexp(  ik.2D* l2D)  *xa(-2D  )d-2D 
Jexp( ik2D-£2D) *AX( k2D)dk2D 


Jexp  (  ^-2d’-2d  ^  *AA^ -2D^-2D 


Jexp  (ik2D*I2D)*xx(k2D)dk2D 


<^*2D)X  ^-2D  +  —20^ 


The  reflected  field  spectrum  above  exhibits  the  same  features 

as  those  of  ocean  surface  roughness  scattering^ ,  ocean  bottom 
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roughness  scattering  '  ,  and  under-ice  roughness  scattering 

The  reflected  field  spectrum  has  two  components.  The  specular 

component  contains  a  factor  of  Dirac  Delta  function  which 

indicates  a  concentrated  energy  in  the  specular  direction.  The 

specular  component  depends  on  the  entire  wave  number  domain  of 

spectrum  associated  with  the  scattering  physical  quantities  such 

as  surface  roughness,  fluctuating  sound  velocity,  etc.  The  rest 

is  the  off-specular  component  which  depends  only  on  a  particular 

(selected)  wave  number  defined  by  the  incident  wave  and  a 

selected  (reflected)  wave  number.  The  specular  component  can  be 

utilized  to  estimate  (forward)  reflection  loss  which  will  be 

reported  in  the  future.  The  off-specular  component  can-be 

utilized  to  estimate  scattering  strength.  The  back-scattering 

strength  in  particular  is  addressed  here.  The  usual  scattering 

2  2 

strength  m"  is  given  by  k  IMkr  2d)  where  vr  is  the  direction 

cosine  of  a  reflected  wave  with  respect  to  the  vertical  axis. 

Accordingly  m"  written  in  terms  of  direction  cosines  (k^z  *  ky, 

k  -  k\>  )  is  given  by  the  following  expression.  (r,v  )  are 
rz  r  it 

direction  cosines  of  (incident,  reflected)  waves  with  respect  to 
the  vertical  axis. 

«'  -4  k^S2'W!ii,2D-!ir,2D>  -  2i46Fxk2/[n<2(r+vr)2+62)- 
•  (k2(Y-vr)2  +  32)  ]  ^XA(ki/2D  -  krf2D)  +  ^pFxk2/[(k2(Y+vr)2  +  e2). 


•(k2w-vr>2*e2u}2*xx(kii2D  -  ktj2Dn 


(2-6) 

The  back-scattering  strength  mB"  (10  log^Q  ms"  in  level)  is 
thus  obtained  by  setting  kr  2d  “  2D  *a*so  ym  vr*  in  t*le  al30ve 

equation.  It  results  in  the  following  expression: 


V  -  4  k4r4|»AAi2iii.2D|-2<4eFXk2'/>|4k2l'2  +  62,B2^V2^,2D» 

+  {  4(3Fxk2/I  (  4k2Y2+e2  )  e2 1  ^2*xx(  2k.  ,  2D)  J  . 


(2-7) 


In  order  to  estimate  backscatter ing  strength,  spectra 

and  P  should  be  estimated  first.  An  example  is 

AA  XA  XX 

worked  out  in  the  next  section  for  demonstrating  its  use. 


3 .  An  Example  and  its  Physical  implications 


The  following  predictions  are  based  on  plausible  assumptions 
of  x  ■  ■  h  and  -  *AX  -  0  .  Physically,  the  wave  breaking 

roughness  elements  are  assumed  in  the  first  assumption  to  serve 
also  as  surface  roughnesses  that  scattered  sound.  In  the  second 
assumption,  r>  and  are  assumed  not  correlated.  This  may  be 
partly  justified  by  the  fact  that  plume  penetration  takes  time 
and  the  local  volume  scattering  bubbles  may  not  have  been 
generated  by  the  local  breaking  roughnesses.  Then  equation  (2-7) 
is  given  by: 


mB"  -  4  k4Y4[i+^eFxk2/n4k2Y2+(B2)e2)  ^2]*llA(2ki  2D) 


(3-1) 


where 

Fx  -  Fnl  -  2na/(«R/c0)2(g/30aCH)(2.9  x  106) 

-  2na(10~6)/f2n  3.25xl06  /1+O.izl2  •  (9.8/30x0.0185)  (2.9xl06) 
c0  a 

*  1.616  X  10-6  a3 

a  -  bubble  radius  in  /jm 

1  +  .lz  -  1.075  average  over  1.5  m  depth 
(3  -  1/1  z  (1/m) 

k  -  0 . 0041867f  (1/m) 


|2k.  2D  I  “  2k/l^V2  -  0.0083734  f  cos9g 

0g  -  grazing  angle 

*AA  (kB)  "  8-lXlO"3/(4RkB4)  exp  [ -0 . 74g2/U4kg2 } 

14 

-  Pierson  -  Moskowitz  Spectrum 

r  -  cose.-  sin0„ 

1  9 

9i  -  incident  angle. 

Substitution  of  above  numerical  numbers  into  equation  (3-1) 
results  in  the  following  equation: 


m  "  -  ( 1+4 1.133  x  10  10f2(al  )3/(l+7.01xl0-5f2  sin©  1  2 ) \2 } • 

O  Z  y  z 

•1.61xlO"4tan40gexp  [ -1 . 01xl06/f 2cos29gU4 ) ] 

-  [1  +  Volumetric  Factor]  (Surface  Scattering  Effect) 

(3-2) 


Accordingly,  the  back-scattering  strength  has  two  independent 
components.  One  component  estimates  the  surface  scattering 
effect  and  has  the  exact  form  utilized  by  Ogden/Erskine*^ .  The 
other  component  estimates  the  volumetric  scattering  effect  and  is 
given  in  this  example  in  terms  of  the  surface  scattering  effect 
multiplied  by  a  volumetric  factor.  This  volumetric  factor  is  the 
sum  of  the  product  of  plume  penetration  effect  and  F  as 
originally  given  in  the  form  of  equation  (2-3)  and  specifically 
for  this  example  it  results  in  equation  (3-1).  The  parameter 


1  -(1/a)  quantifies  the  extent  of  plume  penetration.  The 
parameter  "a"  stress  the  bubble  size  dependence  on  relating 
($c/Cq)  to  other  physical  parameters.  There  is  an  indication 
that  1  depends  on  the  wind  speed.  The  dependence  of  bubble  size 
on  wind  speed  is  less  known.  For  instance,  the  bubble  size 
spectrum  was  concluded  by  Wu*^  to  be  approximately  invariant  with 
either  the  depth  or  wind  speed. 

Before  numerical  computation,  it  will  be  helpful  (1)  to 
review  the  validity  of  (2)  to  investigate  the  qualitative 
behavior  of  equation  (3-2). 

Major  assumptions  made  were  the  applicability  of  perturbation 
method  for  both  volumetric  and  surface  scattering.  For 
volumetric  scattering,  it  amounts  to  assuming  the  validity  of 
single  volumetric  scattering.  It  has  been  estimated*  that  the 
assumption  requires  acoustic  frequencies  to  be  not  higher  than 
1000  Hz.  For  surface  scattering,  the  perturbation  approximation 
requires  maximum  wind  speed  to  be  limited  to  about  13  knot  and  30 
knot  for  1000  Hz  and  100  Hz,  respectively.  Therefore,  the 
following  predictions  for  frequencies  200  Hz  and  lower  are  more 
reliable  at  high  wind  speeds.  Though  the  following  predictions 

5 

are  better  than  those  of  previous  investigators,  McDonald  and 

Henyey®,  who  did  not  include  surface  scattering  effect  and 

utilized  far-field  approximation,  predictions  at  high  frequencies 

around  1000  Hz  for  wind  speeds  exceeding  13  knot  are  considered 

preliminary.  In  the  future,  it  may  be  somewhat  generalized  by 

the  composite-roughness  scattering  concept  (see  e.g.  McDaniel  et 
16  17 

al  ,  Jackson  et  al  ).  Fortunately,  where  the  surface 


scattering  perturbation  is  questionable,  the  volumetric 
scattering  dominates  and  volumetric  perturbation  is  valid. 

The  qualitative  behavior  of  equation  (3-2),  depends 

critically  on  the  magnitude  of  the  volumetric  factor.  It  can  be 

easily  seen  that  the  volumetric  factor  is  much  less  than  one  for 

small  values  of  (f,l  ,a)  where  k  -  2nf.  This  means  that  the 

z 

back-scattering  strength  is  dominated  by  the  surface  scattering 

phenomena  when  frequency  or  wind  speed  is  low.  This  qualitative 

prediction  is  consistent  with  the  experimental  data  summary 

presented  by  Ogden/Erskine1^ .  It  can  also  be  shown  that  the 

volumetric  scattering  effect  dominates  when  1  or  wind  speed  is 

sufficiently  high.  There  is  an  interesting  coincident 

observation.  In  summarizing  the  experimental  data  on 

characteristic  depth,  1  here,  Wu10  concluded  that  1  stayed  at  a 

z  z 

low  constant  value  of  about  0.4m  up  to  about  14  knot  wind  speed 
and  then  suddenly  increased  linearly  with  wind  speed  to  a  1 

z 

value  of  about  lm  at  23  knot.  Ogden/Erskine1^  also  concluded 
that  the  back-scattering  strength  increased  suddenly  beyond  this 
critical  wind  speed  for  higher  frequency  acoustic  waves  -  the 
increase  was  presumed  to  be  caused  by  the  volumetric  scattering 
effect.  These  two  observations  can  be  linked  together  by 
equation  (3-2).  Equation  (3-2)  also  predicts  less  increase  in 
the  volumetric  scattering  effect  for  a  higher  grazing  angle. 

This  prediction  is  also  consistent  with  the  back-scattering  data 
by  comparing  data  for  10°  and  30°  grazing  angles.  Quantitative 
predictions  follow. 


Special  phenomena  of  interest  are  (1)  the  critical  wind  speed 


beyond  which  the  volumetric  scattering  dominates  and  (2)  the 
excess  dB  level  due  to  the  volumetric  scattering.  Following 
estimates  assume  a-60  /jm. 


At  critical  wind  speed,  the  volumetric  factor  becomes  an 

order  of  one.  For  a  given  frequency  and  a  grazing  angle,  a 

critical  lz  can  be  obtained.  Then  the  critical  wind  speed  can  be 

obtained  from  the  expression  describing  1_  dependence  on  wind 

z 

speed.  In  general,  1_  (in  m)  increases  with  wind  speed  (U  in 

z 

ms-* ) .  Wu*^1  proposed  the  following  two  relationships. 

1_  -  0.4m  when  U<7ms-*(«14  knot) 

z 

1  ■  0.4+0.12(U-7)m  when  U>7ms-* 

z 

(3-3) 


and 

1  »  0.85  when  U<7ms-* 

z 

1  -  0 .85+0 . 12 ( U-7 )m  when  U>7ms-1 

z 

(3-4) 


Equation  (3-3)  is  based  on  the  bubble  population  data  of 
18  19 

Kolovayev  ,  Johnson  and  Cooke  and  on  the  acoustic  cross 

20 

section  data  of  Thorpe  .  Equation  (3-4)  is  based  on  the 

21 

acoustic  cross  section  data  of  Crawford  and  Farmer  .  The  wind 

independency  of  1  up  to  7ms-*  wind  velocity  is  a  curious 

z 


physical  phenomenon  that  should  be  investigated  further  in  the 
future.  Up  to  now,  the  temptation  is  to  exclusively  conjecture 
the  phenomenon  as  being  due  to  biological  activities.  There  may 
be  other  explanations  that  are  dynamic  in  nature.  The  estimated 
critical  wind  speeds  for  different  frequencies  and  grazing  angles 
are  tabulated  in  Table  1  and  depicted  in  Figure  2  and  Figure  3. 

The  excess  dB  level  due  to  the  volumetric  scattering  is 
estimated  by  10  log^  [ l+( volumetric  factor)]  and  added  to  the 
surface  scattering  level  10  log^g  [surface  scattering  effect]  to 
obtain  the  total  backscattering  level.  The  predicted  excess  dB 
levels  are  tabulated  in  Table  1.  The  predicted  total 
backscattering  levels  are  depicted  in  figures  2  through  5. 

For  50Hz  and  either  10°  or  30°  grazing  angle,  the  estimated 

critical  wind  speeds  utilizing  either  equation  (3-3),  denoted  by 

Wu,  or  equation  (3-4),  denoted  by  C/F,  are  higher  than  30  knots. 

This  means  that  the  surface  scattering  effect  dominates  and  the 

excess  dB  level  due  to  the  volumetric  scattering  is  negligible 

(zero  in  Table  1)  even  beyond  wind  speed  of  30  knots.  These 

predictions  are  consistent  with  data  summary  of  Ogden  and 
15 

Erskine  (O/E)  as  shown  in  Figure  2  and  Figure  3  in  which  no 
critical  wind  speed  was  observed  for  50  Hz  case  up  to  wind  speed 
of  30  knot.  The  50  Hz  total  backscattering  levels  at  30  knot 
wind  speed  were  estimated  and  depicted  on  Figure  5.  Since  the 
backscattering  levels  are  dominated  by  the  surface  scattering 
effect,  the  term  that  requires  the  use  of  either  Equation  (3-3) 
or  Equation  (3-4)  is  negligible  and  the  estimated  levels  are  both 
identical  to  those  of  Ogden  and  Erskine  shown  in  Figure  5.  The 


same  conclusion  can  be  made  on  the  100  Hz  backscattering  levels 
for  a  given  wind  speed  of  5  knot  as  shown  in  Figure  4. 

For  200  Hz  and  either  10°  or  30°  grazing  angle,  the  estimated 

critical  wind  speeds  and  backscattering  levels  are  closer  to  O/E 

data,  if  equation  (3-4)  or  C/F  is  utilized  as  in  Table  1,  Figure 

2,  and  Figure  3.  The  results  stress  the  importance  of  1  effect 

z 

on  both  critical  wind  speeds  and  the  backscattering  levels. 

Though  the  comparison  of  predictions  and  experimental  data  for 
30°  grazing  angle  case  is  not  as  good  as  that  of  10°  case,  the 
trend  of  lower  excess  dB  due  to  volumetric  scattering  is 
successfully  predicted. 

The  importance  of  lz  effect  on  backscattering  is  further 

demonstrated  by  the  case  of  1000  Hz.  The  use  of  lower  1  values 

z 

given  by  Wu,  predicted  a  critical  wind  speed  of  O/E  for  10° 
grazing  angle  (Figure  1)  but  much  higher  value  for  30°  grazing 
angle  (Figure  3).  The  predicted  excess  dB  levels  above  these 
critical  wind  speeds  are  lower  than  those  observed  by  O/E  in 
general.  However,  the  use  of  higher  1  values  given  by  C/F 
predicted  no  critical  wind  speed  for  either  10°  or  30°  grazing 
angle  case  -  i.e.  the  volumetric  scattering  dominance  is 

predicted  at  wind  speed  as  low  as  6  knots.  In  Figure  2,  the  data 

22  2  3 

of  Christian  and  Tattersall  (C/T),  Chapman  and  Harris  /Chapman 
24  25 

and  Scott  (CH/CS),  and  Chester  ( CST7 )  are  presented.  Though 
these  data  are  for  7°  grazing  angle,  they  clearly  indicate  no 
critical  wind  speed  above  6  knots.  The  trend  is  similar  to  that 
predicted  by  utilizing  higher  1  values  of  C/F.  If  1  data  of 
C/F  are  indeed  mixed  with  bubble  population  of  biological 


1  0 

activities  (beside  wind  wave  generated)  as  speculated  by  Wu  , 

then  the  wave  mixing  tends  to  produce  more  scattered  1 data  as 

2 1 

indicated  by  the  original  data  of  Crawford  and  Farmer  .  This 

implies  that  the  spread  of  backscatter ing  levels  in  the  region  of 

Figure  2  between  those  predicted  by  C/F(1000)  and  Wu(1000)  is 

quite  possible.  As  expected,  the  predicted  grazing  angle 

dependence  at  low  wind  speed  for  1000  Hz  (Figure  4)  compares 

quite  well  with  that  of  O/E  data  if  lower  1  values  of  Wu  are 

z 

used  (when  C/F  1  values  are  used,  levels  are  closer  to  C/T 

Z 

data).  At  high  wind  speed  of  30  knots  (Figure  5),  the  predicted 
grazing  angle  dependence  compares  well  with  O/E  data  only  at  low 
grazing  angles  and  if  higher  1  data  of  C/F  is  utilized. 

*  Z 


4.  Conclusions  and  Recommendations 


The  preliminary  analysis  based  on  selecting  roughness  element 

height  as  a  choice  physical  parameter,  succeeded  in  predicting 

observed  features  of  backscattering  strength  levels,  such  as  a 

critical  wind  speed  beyond  which  volumetric  scattering  dominated. 

It  was  especially  interesting  that  this  critical  wind  speed 

corresponded  to  the  wind  speed  beyond  which  1  suddenly  increased 

linearly  with  the  wind  speed.  It  was  also  found  that  wind  speeds 

alone  could  not  uniquely  define  1  and  that  quite  different 

z 

backscattering  strength  levels  were  predicted  for  different 
experimental  data  set  taken  at  the  same  apparent  wind  speed.  If 
wind  speed  were  to  be  used,  one  should  use  additional  parameters 
such  as  wind  fetches  and  durations.  Otherwise  one  may  use  local 
physical  parameters  such  as  wave  heights  or  friction  velocities 

O 

suggested  by  Wu  . 

Similar  analyses  based  on  choices  of  other  physical 
parameters  should  be  performed  to  relate  scattering  strength  with 
those  parameters,  e.g.  Nq,  W,  and  Myg*  These  analyses  require 
information  on  their  cross-correlations  or  their  wave  number 
spectra  of  these  chosen  physical  parameters.  It  will  be  better 
if  those  physical  parameters  are  measured  simultaneously  with 
backscattering  data.  CST-7  was  such  an  experiment  and  its  data 
should  be  utilized.  Then  the  prediction  of  reflection  loss 
should  be  performed.  If  successful,  the  same  mathematical  method 
can  be  utilized  to  evaluate  ocean  bottom  scattering  when 
volumetric  scattering  is  important. 
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TABLE  1  Backscattering  Phenomena  and  e-folding  bubble 
Plume  Depth  (lz): 

Vu  =  (Kolovayev,  Johnson/Cooke,  Thorpe) 

C/F  =  Cravford/Farmer ,  0/E  =  Ogden/Erskine, 
C/T  =  Christian/Tattersall 
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L.  Dillman,  R.  Nielsen,  R.  Saenger,  M.  Sundvik,  J.  Tattersall, 

M.  Vaccaro 

3111  D.  Klingbeil,  S.  Capizzano,  J.  Turner,  T.  Baus 

312  M.  Gerber 

3121  R.  Dullea,  N.  Fisch,  W.  Powers,  D.  Browning,  M.  Rosario 

3122  E.  Jensen,  J.  Bishop,  H.  Weinberg 
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33  C.  Nawrocki 

33A  P.  Berstein,  A.  Goodman,  P.  Koenigs,  J.  Monti,  R.  Malone, 

L.  Petitpas 

3314  R.  Barton,  R.  Tremblay 

332  M.  Ricciuti 

3321  P.  Seaman,  E.  Kuo  (5) 

61  C.  Batta,  D.  Williams 

NEWPORT 

Code 

0262  Library 
10  K.  Lima 

8212  G.  Bishop 


External:  86 
Internal:  59 
Total:  175 


